Equilibration and Dynamic Phase Transitions of a Driven Vortex Lattice 
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We report on the observation of two types of current driven transitions in metastable vortex lattices. 
The metastable states, which are missed in usual slow transport measurements, are detected with 
a fast transport technique in the vortex lattice of undoped 2H-NbSe2. The transitions are seen by 
following the evolution of these states when driven by a current. At low currents we observe an 
equilibration transition from a metastable to a stable state, followed by a dynamic crystallization 
transition at high currents. 



PACS numbers; 74.60. Ge 74.60. Jg 74.60Ec 

Vortices in type II superconductors are easily trapped 
in long lived metastable states created by the random 
pinning potential A driving current assists in vor- 
tex de-trapping and in finding the stable states |§-D, 
in essence assuming the role of temperature in ordi- 
nary phase transitions. This analogy was pointed out 
by Koshelev-Vinokur (KV) who showed that when the 
vortices are set in motion by a current /, the random 
potential appears as a temporally fluctuating Langevin 
force in the moving frame ^ . The random potential can 
then be replaced by a "shaking temperature" Tg cx l/I, 
which leads to a simpler problem of a pure system at 
an effective temperature T^ff = T + Tg. Thus, if the 
pure system crystallizes at a temperature T,„, the ac- 
tual system will crystallize when driven by a current 
/( oc {Tm — T)~^. This KV transition separates a disor- 
dered state at low currents from an ordered one at high 
currents. Further work ||7[-p^ predicted that the transi- 
tion is preceded by a regime of plastic flow followed by 
smectic ordering. Small angle neutron scattering (SANS) 
on the field-cooled (FC) vortex lattice in undoped 211- 
NbSe2 [ p^ provided evidence for current induced order- 
ing. This data are generally accepted as direct evidence 
for the KV transition [p| Jll|Jl^ despite the fact that for 
usual slow transport measurements under the same ex- 
perimental conditions the transition is absent ||l^,|l^,|l6| . 
This discrepancy could be an indication of metastability 
in the FC vortex state, in which case the observed order- 
ing in the SANS data would be an equilibration rather 
than the KV transition. Thus far however metastability 
in the vortex lattice of undoped 2H-NbSe2 was thought 
to be absent JlSf . 

In this letter we present results of fast transport mea- 
surements that demonstrate the existence of metastable 
states in undoped 2H-NbSe2 . The experiments probe the 
dynamic transitions with emphasis on distinguishing be- 
tween equilibration and the KV-transition. We show that 
the current induced transition in the FC state is in fact 
an equilibration from a disordered metastable state to an 
ordered stable state and not the KV transition. We fur- 



ther find that in the peak effect region (a region of the 
phase diagram, just below Hc2 where the critical current 
increases with field and temperature ||]) the vortex lat- 
tice can undergo two current induced transitions. The 
first is an equilibration transition observed at low cur- 
rents which is then followed by the KV crystallization at 
much higher currents. 

The sample was an undoped single crystal 2H-NbSe2 
platelet of dimensions 1.5x0. 65x0. 025mm'^. Its critical 
temperature Tc was 7. IK , the transition width 80 mK 
and the normal resistance near Tc was 21xl0~'^il. A 
four probe measurement with low resistance Agin sol- 
der contacts was used to monitor vortex response. The 
response to fast current ramps - 200 A/s- was detected 
with a fast (2/l(s response time) amplifier while the slow 
ramp -5 x 10~^A/s - measurements were obtained with 
a Keithley 181 nanovoltmeter. The magnetic field was 
kept along the c axis of the sample and the current was 
in the a ~ b plane. The zero-field-cooled (ZFC) and the 
FC vortex lattices were prepared by applying the mag- 
netic field after or before cooling the sample through Tc 
respectively in the absence of applied current. The de- 
gree of order of the vortex lattice was inferred from the 
critical current, Ic, (defined by a 5/iV response criterion) 
by using the Larkin-Ovchinnikov model [0, according 
to which Ic grows with the degree of disorder. 

In Fig.l we compare the current-voltage {I — V) curves 
of the FC and ZFC lattices. When probed with slow cur- 
rent ramps the response of the two states is identical 
and no evidence of metastability is observed. This is in 
contrast to results in Fe doped 2H-NbSe2 where strong 
metastability and hysteresis were found in slow transport 
measurements [Q. Here the temperature dependence of 
the critical currents (inset of Fig. 1(a)) is the same for 
both states, and exhibits a pronounced peak effect just 
below Tc- But in spite of the identical response in slow 
measurements, the initial vortex states prepared by ZFC 
and by FC are not identical. The difference between the 
two states becomes evident only in fast measurements. 
This is illustrated in Fig.l (b) where significant hystere- 
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FIG. 1. (a) 7 - 1/ curves obtained with slow (5xlO"^A/s) 
current ramps. The inset shows the temperature dependence 
of the critical current for ZFC and FC vortex lattices obtained 
with slow measurements, (b) Fast I — V curves (200A/s) for 
an FC vortex lattice. The inset shows the time evolution of 
the response of an FC vortex lattice to a current step. 



sis is seen when comparing the fast I ~ V curve for the 
pristine FC state, measured on the first ramp-up of the 
current, with that recorded when the current is ramped 
down. The critical current of the pristine FC state is 
almost twice that of the annealed state recorded on the 
down ramp, indicating that it is more disordered. By 
contrast the pristine ZFC state starts out with a low 
critical current, is unchanged by current cycling and its 
I — V curves are independent of ramping speeds. We con- 
clude that the ZFC state is ordered and stable. For the 
FC state, the J — V curves obtained after the first ramp 
up (including the ramp down and all subsequent ramps) 
are identical to those of the ZFC state indicating that 
the FC lattice reorders during the first ramp. It follows 
that the FC lattice is initially in a disordered metastable 
state which reorganizes under the influence of a current 
into a stable ordered configuration. 

The current driven organization of the FC lattice is 
seen directly in the inset of Fig. 1(b) through the evolu- 
tion of the response to a current step. After an initial 
waiting time the response starts growing from zero as 
the vortices order into a state with lower critical current, 
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FIG. 2. (a) Critical currents for initial and annealed vortex 
lattices prepared by ZFC and FC processes following a 100s 
drive with various currents, I. The open circles are the I — V 
curve of an FC vortex lattice obtained in a slow (5xlO-^A/s) 
measurement, (b) Temperature dependence of /* and the 
critical currents of the FC lattice for fcist and slow measure- 
ments. 



and saturates to a value that depends on the amplitude 
of the applied current. In order to study the current 
dependence of the reorganization the pristine FC state 
was driven with long (lOOsec) current steps and then 
quenched by suddenly removing the current. When the 
current is removed the evolution of the vortex state is 
arrested instantaneously and then resumes from the ex- 
act same state when the current is turned back on. This 
is clearly seen in the inset of Fig. 1(b). Thus, when the 
quenched state is probed by recording the 1 — V curve 
with a fast current ramp the critical current reflects the 
degree of order of the state at the moment of quenching. 
In Fig. 2 (a) we plot the critical current of the quenched 
state as a function of current-step amplitude. No evo- 
lution is seen until I* =6. 7mA where we note the onset 
of a sharp transition from disordered state (higher /c) to 
ordered state (lower 1^). The transition is completed at 
/cs =8 mA, at which point the voltage response becomes 
measurable in the l—V curve. Also shown are the results 
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FIG. 3. (a) Evolution of the critical current with driving 
current amplitude for ZFC, FC and annealed vortex lattices in 
the peak effect region following a 100s drive with continuous 
(CC) and pulsed (PC) currents, (b) Current dependence of 
differential resistance illustrating the definition of Ip and It- 
The data saturates at Rbs, the Bardeen-Stephen free flux 
flow value. Also shown is the I — V curve obtained with a 
pulsed current. 



for the same experiment on the pristine ZFC lattice and 
the annealed FC and ZFC lattices. The vortex lattice 
was annealed with a slow (5 x 10~^A/s) cycle of the cur- 
rent between and 50 mA. The absence of a jump in Ic 
indicates that the current does not cause ordering of the 
ZFC and annealed states. For higher amplitude current 
steps, / > Ics the I — V curves of all quenched states are 
identical, regardless of the initial state. The transition at 
/* is the same as that seen in the SANS measurements 
and interpreted as the KV crystallization But the 
fact that it occurs at such low currents makes it an un- 
likely candidate for the KV transition. This is confirmed 
by the temperature dependence of /*, shown in Fig. 2(b). 
Here /* decreases with increasing temperature, in con- 
trast to the predicted increase with temperature for the 
KV-transition. We conclude that /* is not where the KV 
transition occurs but rather the onset of a current driven 
equilibration from the metastable disordered FC state to 
a stable ordered state. 

To further illustrate the difference between the equi- 
libration and the KV transition we repeated the above 
experiment in the peak effect region, where the KV tran- 
sition ||l5| , p^ and a current driven reorganization in the 
the ZFC lattice H were seen with transport measure- 
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FIG. 4. Temperature dependence of the currents sep- 
arating the different dynamic regimes: I < I* - pinned; 
I* < I < Ics - equilibration; Ic < I < Ip - plastic flow; 
Ip < I < It - smectic ordering; I > It - crystalization. The 
dashed an dotted curves show fits to the K-V model discussed 
in the text. 



ments. Since much higher currents are required in order 
to observe the KV transition. Joule heating needs to be 
considered. For example, for a 45 mA (the highset cur- 
rent in Fig. 3) continuous current drive, a temperature 
increase of 25 mK in the sample was detected. The de- 
tails on the determination of Joule heating in these ex- 
periments is presented elsewhere To avoid heating 
at high currents we apply the current step in a sequence 
of short, 10/is, pulses separated by 500/Lts cooling inter- 
vals with no current. For currents below 15mA pulsing 
was not necessary since heating was negligible (less than 
3mK) and the measured Ic was the same for continuous 
and pulsed applications. In the following discussion we 
focus on the data which show no heating. In Fig. 3(a) we 
show the critical current of the quenched vortex lattice 
following the application of 100s current steps of various 
amplitudes. At low currents / < /*=4.5mA no mea- 
surable change occurs in the vortex states. For higher 
currents, I* < I < /cs=8.6mA the current driven organi- 
zation sets in, affecting each state differently. The critical 
current of the FC state drops rapidly, that of the ZFC 
state increases while the annealed lattice curve is almost 
unchanged. At Ics all the curves converge, overlapping at 
higher currents, which indicates that in this regime the 
vortex state is determined by the driving current alone 
and is independent of the initial preparation. The critical 
current decreases with increasing current step amplitude 
saturating at a value which is lower than that of the sta- 
tionary vortex lattices, as expected for the KV transition. 
In order to identify the KV crystallization we follow 



the usual approach in transport measurements 18 
focusing on the shape of the differential resistance curve 
shown in Fig. 3(b). The differential resistance is vanishly 
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small for low currents and becomes finite only above Ics 
where metastability has disappeared. As the current is 
increased further the curve rises to a maximum value of 
21.7 mfi reached at Ip = 20.4mA and then drops down 
saturating at the Bardeen-Stephen free flux flow value 
Rbs= ISmf] = RnH/Hc2 at It =28. 8mA. Its shape can 
be interpreted according to numerical simulations on the 
motional organization of a vortex lattice 
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low velocities the motion is plastic leading to a defective 
lattice. The defect density increases with velocity result- 
ing in a corresponding increase in differential resistance 
which peaks at Ip. At this point the vortices move in an 
array of almost periodic channels forming a smectic state 
characterized by transverse order alone. As the velocity 
is further increased the vortices order inside the channels 
eventually crystallizing at It- The temperature depen- 
dence of It, shown in Fig. 4, is consistent with the KV 
prediction for the crystallization with r,„ = 4.37±0.01K. 
In the same figure we also show the curves for J*, Ics, Ip 
and It which separate the different dynamic regimes of 
the moving vortex system as described in the caption. 

In both sets of measurements presented here the FC 
state is initially metastable and can reorganize into a 
stable state when driven by a current. The reorganiza- 
tion sets in at a current /* which is significantly lower 
than the critical current of the pristine FC state. Below 
the peak region, the FC state undergoes an equilibration 
transition at I* which leads to an ordered stable state 
indistinguishable from the ZFC state. By contrast the 
ZFC state is unaffected by the current. For the data in 
the peak effect region, vortex reorganization sets in at /* 
for both the FC and ZFC states, one becoming more or- 
dered the other more disordered until, at Ics, they merge 
into a single state and loose the distinction of their ini- 
tial preparation. Thus, by analogy with the equilibration 
transition seen below the peak region, we conclude that 
in this part of the phase diagram the initial FC and ZFC 
states are both metastable and /* is the onset of a current 
driven equilibration transition to a mixed stable state. 

Several aspects of this data can be understood in the 
context of recent Hall probe measurements on Fe doped 
samples |^ which have shown that in the peak effect 
region a stable disordered phase can be injected through 
a surface barrier at the sample edges. This mechanism 
can readily account for the disordering of the ZFC vortex 
lattice as indicated in Fig. 3 and in Ref. ||]. In fact it 
is possible that /* is associated with a surface barrier 
PH , so that bulk vortices are not subject to a Lorentz 
force until /* is exceeded. Other features of the data can 
be described in terms of the phenomenological parameter 
for the distance that the disordered phase penetrates into 
the sample, (T, H, I) . For example — below the 
peak where the ordered state is stable , but it is finite 
at low velocities in the peak region where the two phases 
coexist. At high velocities, where the random potential is 
averaged out in the moving vortex frame, Lj. — Q which 



leads to a stable ordered state. 

In summary, the experiments described here demon- 
strate that a metastable vortex lattice can undergo two 
types of current driven transitions. Below the peak effect 
region we observed a current driven equilibration transi- 
tion from a metastable disordered state to a stable or- 
dered one. In the peak effect region we detected both 
the equilibration and the KV transitions. 
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